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q ' Abstract. We have analyzed the intermediate term periodicities in soft X-ray flare 

, index (FIsxr) during solar cycles 21, 22 and 23. Power spectral analysis of daily 

FIsxr reveals a significant period of 161 days in cycle 21 which is absent during 
cycle 22 and 23. We have found that in cycle 22 periodicities of 74 and 83 days are 
. in operation. A 123 day periodicity has been found to be statistically significant 

during the part of the current solar cycle 23. The existence of these periodicities has 
been discussed in the light of earlier results. 

<n : 

> ■ 

' 1. Introduction 

t> : 

. In the solar cycle 21, an intermediate term (or midrange) periodicity 

of 154 days was discovered by Rieger et al. (1984) in gamma ray and 
soft X-ray flare data. Since then a number of studies have been made to 
look for the evidence of near 155 days periodicity in solar activity. With 

' solar flares and flare related data it has been detected by Bogart and 

9^ Bai (1985), Ichimoto et al. (1985), Ozgiic and Atag (1989), Droge et al. 

(1990), Bai and Oliver (1990), and Kile and Oliver (1991). The presence 

■ of this periodicity was also analyzed in several other solar activity fea- 

tures, such as, sunspot area (Lean, 1990; Carbonell and Ballester, 1990, 
1992), sunspot blocking function (Lean and Brueckner, 1989), 10.7 cm 
•*~J ■ radio flux (Lean and Brueckner, 1989), and Photospheric magnetic flux 

rN \ (Ballester, Oliver, and Carbonell, 2002). On the other hand no evidence 

for the presence of a near 155 days periodicity has been found in any 
solar activity indicator during solar cycle 22 and part of cycle 23 (Kile 
and Oliver, 1991; Ozgiic and Atag, 1994; Oliver and Ballester, 1995; 
Ballester, Oliver, and Carbonell, 2002; Bai, 2003) instead some other 
prominent periodicities were detected in these studies. 

In the present paper, we have made an attempt to investigate the 
existence of intermediate term periodicities in the soft X-ray flare in- 
dex (FIsxr), which is based on the continuous record of soft X-ray 
(SXR) flares observed by GOES during solar cycles 21, 22, and 23. 
The FIsxr is calculated by weighing peak intensity of SXR flares of 
different classes (B to X) in units of 10 -6 Wm~ 2 (see Joshi and Joshi, 
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Figure 1. 27-day moving average of daily FIsxr for the period 1976—2004. The 
dashed line shows the smooth spline curve. 



2004 and references therein). In this manner the daily values of FIsxr 
represent the daily SXR flare activity and are suitable for a short term 
solar activity analysis. To compute the FIsxr, the SXR flare data for 
the time span of 01 June 1976 to 31 May 2004 has been downloaded 
from NGDC's anonymous ftp server: ftp://ftp.ng dc.no aa.gov/STP/ 
SOLAR_DATA/SOLAR_FLARES/XRAY_FLARES. In Figure 1 the 
time variation of 27-day moving average of daily FIsxr for the period 
1976-2004 is presented. 

To perform the study, we have split the data in three parts corre- 
sponding to the period of solar cycle 21 (June 1976 to August 1986) 
with 3717 days data, cycle 22 (September 1986 to April 1996) with 
3524 days data and cycle 23 (May 1996 to May 2004) with 2953 days 
data. Cycle 23 is still not complete and the data for this cycle covers 
the ascending phase, the maximum and part of the descending phase of 
the solar cycle. Also, for cycle 23 we have performed the power spectra 
for three different time intervals corresponding to the periods of time 
studied by Bai (2003). 



2. Periodogram Analysis 

We have Used the Lomb-Scargle periodogram method (Lomb, 1976; 
Scargle, 1982) modified by Home and Baliunas (1986) to perform the 
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Figure 2. (a) Normalized power spectrum of the daily FIsxr time series of solar 
cycle 21 for the frequency interval of 231.5—57.8 nHz (50—200 days). Horizontal 
dashed lines indicate FAP levels of 10 %, 30 %, and 50 % (from top to bottom) (b) 
Scargle power distribution corresponding to (a) . The vertical axis is the number of 
frequencies for which power exceeds x. The straight line is the fit to the points for 
lower values of power. 



power spectra of FIsxr time series. The power spectra for the three 
cycles 21, 22 and, 23 has been computed for the frequency range of 
231.5—57.8 nHz which corresponds to the period interval of 50—200 
days. Figure 2(a) shows the normalized power spectra of daily FIsxr 
values for solar cycle 21. The FIsxr is not independent. Therefore the 
probability P of the power density at a given frequency being greater 
than K by chance is given by 

P(z > K) = exp(-K/k), (1) 
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Figure 3. Same as Figure 1, but for solar cycle 22. 

where the normalization factor k, which is due to event correlation, can 
be determined empirically (Bai and Oliver, 1990). 

Figure 2(b) shows the distribution of the power values corresponding 
to the normalized power spectra shown in Figure 2(a). The vertical 
axis shows the cumulative number of frequencies for which the power 
exceeds a certain value. For example, in the power spectrum for solar 
cycle 21 shown in Figure 2(a) we have 174 frequencies. For all these 
frequencies the power exceeds zero; thus, we have a point at (X=0, 
Y=174). At only one frequency (72.0 nHz, which is equivalent to 160.75 
days) the power was 27.42, its maximum value. For lower values of 
power, the distribution can be well fitted by the equation y = 174 
exp(-x/4.34), as expected from equation (1). Thus, we normalize the 
power spectrum by dividing the powers by 4.34 to obtain Figure 2(a). 
For other cases we use the same procedure for normalization. 
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Once the power spectrum has been normalized properly, we can use 
the 'False alarm probability' (FAP) formula to estimate the statisti- 
cal significance of a peak in the power spectrum. It is given by the 
expression 

F = 1 — [1 — exp(-Z m )] N , (2) 

where Z m is the height of the peak in the normalized power spectrum 
and N is the number of independent frequencies. The interpretation 
of F is as follows: if we have a discrete power spectrum giving the 
power at each of N independent frequencies for a set of random data, 
F indicates the probability that the power at one or more of these 
frequencies will exceed Z m by chance. Fourier components calculated 
at frequencies at intervals of independent fourier spacing, A/j^ a =r -1 , 
where r is the time span of the data, are totally independent (Scargle, 
1982). For example, in the power spectrum shown in Figure 2(a), we 
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Figure 5. Same as Figure 1, but for the interval from May 1996 to October 2002 of 
the solar cycle 23. 



have r=3717 days and A/j/ s =3.1 nHz. Thus there are 56 independent 
frequencies in 231.5—57.8 nHz interval. We have oversampled to ob- 
tain this power spectrum in which the height of the peak at 160.75 
days is 6.3. The oversampling tends to estimate more accurately the 
peak value. Therefore, if we substitute Z m =6.3 and N=174 (since we 
searched 174 frequencies with 1.0 nHz intervals) in equation (2) we 
get F=0.27, i.e., the probability to obtain such a high peak at 72.0 
nHz (160.75 days) by chance is about 27%. The same analysis has 
been applied to the other power spectra also shown in Figure (3), (4), 
and (5). In the power spectra for solar cycle 22 (Figure 3) we have 
found two important peaks at 157.0 nHz (73.72 days) with FAP=19 
% and 140.4 nHz (82.67 days) with FAP=32.9 %. No significant peak 
appears in the power spectra for solar cycle 23 (Figure 4) when the 
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data from the beginning of cycle 23 to 31 May 2004 was analyzed. We 
have also performed the power spectra for three different time intervals 
(May 1996 to October 2002, September 1999 to June 2001 and October 
1999 to April 2001) studied by Bai (2003). We could not find any 
statistically significant peak in the second and third duration. However 
in the first time interval (i.e., May 1996 to October 2002) a significant 
peak appears at 94.0 nHz (123.13 days) with FAP=11% (Figure 5). 

To be sure that peaks in the different power spectra (Figures 2, 
3, and 5) are not due to aliasing, we have removed sine curves of 
their periods from the original time series. We see that these peaks get 
removed in the power spectra of time series obtained after subtraction. 



3. Summary and Discussion 

We have studied periodicities in FIsxr during solar cycles 21, 22 and 
23 separately. It will be useful to compare and discuss our results in 
the light of periodicities detected in other solar activity indicators. 

We have found that a period of 161 days appeared in FIsxr during 
solar cycle 21, which is absent during cycle 22 and 23. The near 155 
days periodicity, first detected by Rieger et al. (1984) in gamma ray 
and soft X-ray flares during cycle 21, has been found in several other 
studies. During the same solar cycle, it was detected in flare producing 
energetic interplanetary electrons (Droge et al., 1990), proton flares 
(Bai and Oliver, 1990), and microwave flares (Bogart and Bai, 1985; 
Kile and Oliver, 1991). This periodicity was also found in Ha flares 
(Ichimoto et al., 1985), microwave flares (Bogart and Bai, 1985), and 
Ha flare index (Ozgiic and Atac, 1989) when the data of solar cycles 20 
and 21 was considered together. Lean and Brueckner (1989) detected 
a periodicity at 155 days in the sunspot blocking function and the 
10.7 cm radio flux during solar cycles 19, 20 and 21. Carbonell and 
Ballester (1990, 1992) have shown that a periodicity around 150—160 
days seems to be significant during all solar cycles from 16 to 21 in 
sunspot areas. Recently Ballester, Oliver, and Carbonell (2002) in the 
analysis of historical record of photosphere magnetic flux detected a 
significant periodicity near 160 days during solar cycle 21. 

Several mechanism have been put forward to explain the cause of 
near 155 days periodicity. Bai and Sturrock (1991) and Sturrock and 
Bai (1992) proposed that Sun contains a "clock" with a period of 25.8 
days (later modified as 25.5 days) and suggested that periodicity of 154 
days is a subharmonic of that fundamental period. Bai and Sturrock 
(1993) analyzed the longitude distribution of major flares of cycles 
19—22 and interpreted this hypothetical clock as being an obliquely 
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rotating structure (or a wave pattern) rotating with a period of 25.5 
days about an axis titled by 40° with respect to the solar rotation 
axis. On the other hand, another explanation suggests that this peri- 
odicity is associated with the periodic emergence of the magnetic flux 
linked to the regions of strong magnetic field (Carbonell and Ballester, 
1990, 1992; Oliver, Ballester, and Baudin, 1998; Ballester, Oliver, and 
Baudin, 1999; Ballester, Oliver, and Carbonell, 2002). 

Other periodicities detected by us in FIsxr during cycles 22 and 23 
are also well consistent with the earlier work. We detected 74 days and 
83 days periodicities in FIsxr for cycle 22. Bai (1992) analyzed the 
major flare occurrence of solar cycle 22 and found a periodicity of 77 
days which was interpreted as the third subharmonic of the 25.5 days 
period. Ozgiic and Atag (1994) also found a period of 73 days in Ha 
flare index during cycle 22. Oliver and Ballester (1995), in the analysis 
of sunspot areas, found a periodicity around 86 days (close to 83 days 
period in FIsxr) which is statistically significant during part of solar 
cycle 22. For cycle 23 we have carried out analysis in four different 
time intervals, but a significant peak appears only in the duration of 
May 1996 to October 2002 with a period of 123 days. Lou et al. (2003) 
found a periodicity of 122 days in X-ray flares of class > M5.0 during 
the maximum phase of solar cycle 23. Bai (2003) detected a periodicity 
of 129 days in major flares during part of cycle 23 (September 1999 
to June 2001). Here it is interesting to note that 74 and 123 days 
periodicities, found in FIsxr, are close to integral multiple of 25.5 
days. 
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